Dileptons from hot heavy static photons 
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Abstract 

We compute the production rate of lepton pair by static photons at finite temperature at two-loop 
order. We treat the infrared region of the gluon phase space carefully by using a hard thermal loop gluon 
propagator. The result is free of infrared and collinear divergences and exhibits an enhancement which 
produces a result of order ~ e 2 g 3 instead of ~ e 2 g 4 as would be expected from ordinary perturbation 
theory. 



1 Introduction 

The production of photons and leptons at finite temperature is studied as a signal for the formation of a 
quark-gluon plasma in heavy ion collisions ^El El El El El El ■ Photons and leptons are expected to escape 
easily from nuclear matter because of their large mean free path. Computations are done in two steps. The 
first step is to compute the production rate at a given order in perturbative thermal field theory. The second 
step is to substitute this result into a hydrodynamical model that simulates the experimentally realized heavy 
ion collision 01 El El El- In this paper we discuss only the first of these steps. 



The production rate of lepton pairs by a photon of virtuality Q 2 is [5] 



dN 



a 1 

Y2^3 Q2 nB W° ) Im U R (90 , q) ^ 



(1.1) 



dqo d 3 qd 4 x 

Note that there are two coupling constants involved in this equation, or equivalently, two fine structure con- 
stants: e or a characterizes the electromagnetic interaction, and g or a s characterizes the strong interaction. 
The formula above is valid to first order in a and arbitrary order in a s , although the actual calculation will 
involve a perturbative expansion of the polarization tensor. 

In this paper we are interested in contributions to lepton pair production by heavy static photons. We 
take the photon momentum to be: Q M = (<2o>0) with Q 2 = 2> g 2 T 2 . At the Born level, the imaginary 
part of the trace of the polarization tensor is given by : 



lmn fl ( 9o ,0)^ =Q N c q 2 



1 



2 ™F ( y 



(1.2) 



which we denote, for simplicity, ImII Born . N c is the number of colors. 

The rate of static lepton pair production has also been calculated at the two-loop level in the bare 
theory. In physical terms these corrections include the contributions of Compton scattering, quark-antiquark 
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annihilation into a virtual photon, and the 3 
been found to be 151: 



1 process qqG — > 7*. In the large q^/T 2 limit, the result has 



Imn K (go,0)"„ =Imn E 



4m 2 

~9o~ 



2m 2 



g 2 ' ° f 44^ 



(1.3) 



where m 2 



Cf9 2 T 2 /A is the thermal quark mass generated by the one-loop self-energy corrections, with 
Cf = 4/3 the usual QCD structure constant. The temperature independent term is simply the T = 
first order QCD correction to the process e + e _ — ► qq with massless kinematics. The square root threshold 
factor in eq. i|1.3|) is obtained if the thermal self-energy corrections are resummed on the quark lines. The 
interpretation of the temperature dependent corrections is very simple. Starting from massless quarks the 
thermal corrections, in the considered limit, are obtained by simply considering the production in the vacuum 
of quarks with a thermal mass. This result supports the quasi-particle interpretation of thermal effects. In 
the strict perturbative sense, which we adhere to here, the threshold factor can be expanded in to 2 /<7q and 
one observes a cancellation of thermal corrections at 0(g 2 ): the first thermal corrections are then expected 
to appear 0(g 4 ), 

3 q 2 „,„r 4 " 



Imn fi (g ,0)^=Imn 



Born 



i + c F 



4 4tt 2 



4) 



(1.4) 



These calculations are done by integrating over the full fcrmion and gluon phase space including con- 
tributions from both hard and soft modes. However, it is well known that it is inconsistent to use bare 
propagators at soft momentum scales {p ~ gT). For soft momenta, collective modes play an important 
role and hard thermal loop (HTL) resummed propagators must be used ETJIE^E^]. We expect that this 
problem is more serious for gluon modes than for fermion modes. Gluons are bosons and carry a thermal 
factor of the Bose-Einstein form which has the limiting form T/Iq ~ 1/j > 1 when Iq (the gluon energy) is 
~ gT. In contrast, fermions carry a thermal factor of the Fermi-Dirac form which approachs 1/2 in the high 
temperature limit. Thus we expect that collective effects due to soft gluons give rise to an enhancement. 

In this article we will calculate the two loop photon polarization tensor using a HTL gluon propagator 
in order to study the importance of the soft part of the gluon phase space (Fig. pQ). When the resummed 
gluon propagator is used the physics is changed in the following ways: 

[1] transverse and longitudinal modes have different behaviour and must be considered separately; 
[2] in the time-like region gluons become massive quasi-particles; 

[3] in the space-like region the gluon polarization tensor has a non-zero imaginary part which leads to Landau 
damping. This damping is associated with new physical processes which do not appear at the bare two loop 
level and correspond to the scattering of quarks and gluons in the medium. 

One interesting aspect of the study of photon and lepton pair production at finite temperature is the 
treatment of infrared and collinear divergences. At zero temperature there is a general result (the KLN 
theorem E3J E31) which insures the cancellation of these divergences. At finite temperature, the problem 
is complicated by the presence of the collective modes. In all calculations divergences are removed through 
cancellations that are similar to those that occur at zero temperature, however, there is no generalized version 
of the KLN theorem that applies at finite temperature. 

In bare perturbation theory, it has been shown that the corresponding diagrams are finite with respect 
to collinear and infrared divergences ^3 E2J Ell EH EH- The case of a non-static photon has been studied 
in 1201 12] • A calculation using resummed propagators has been done in 22 . The authors use a scalar 
external particle, bare fermion propagators and a HTL gluon propagator. They find that all divergences 
cancel, and that the contribution from longitudinal gluon modes is enhanced relative to the contribution from 
transverse gluon modes by one power of the coupling constant. The transverse mode gives a contribution to 
the polarization tensor of order e 2 0(g 4 ), as is expected from perturbation theory, and the longitudinal mode 
gives a contribution of order e 2 0(g 3 ). In this paper, we perform the corresponding calculation for a virtual 
photon (for which the result is related to the dilepton production rate l|l.lfl which is a physical quantity). 
We find that the cancellation of divergences has the same structure as in the scalar theory, but we obtain an 
enhanced result ~ e 2 g 3 in both the longitudinal and transverse sectors. 

This paper is organised as follows : in Sec. |2]we review some elements of real time finite temperature 
field theory and of the hard thermal loop effective theory. In Sec. we outline the calculation of the 2-loop 
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diagrams. We discuss some preliminary results in Sec.^J In Sec. we perform the final integrals over the 
gluon energy and momentum. Finally, in Sec. we present some conclusions. 

p 




(a) (b) (c) 

Figure 1 : Vertex (a) and self-energy ((b) and (c)) diagrams 



2 Hard thermal loops 

We consider the case of hard fermions and soft gluons. Fermion propagators have momenta of order p 3> gT, 
and are not modified by collective effects. The fermionic propagator is given by 

i 



S(P)=fD(P); D{P) = -s 

P Po - P 



Retarded and advanced functions are obtained by taking the appropriate ie prescription 



D(P) - D R/A (P) = „ 



n2 



We also define the spectral function and the principal part, 

p{P) = D R (P) - D A (P) = 2n c(po) 5{P 2 ) = -- V n S( Po + np) (2.1a) 

V{P) = 2 ( (p + it? -p 2 + (po - ie)* - p 2 ) (2 ' lb) 

which gives, 

D R (P) = iV(P) + \ P {P) (2.2a) 

D A {P)=iV{P)-\p{P) (2.2b) 

The exchanged gluon is soft, and thus the propagator will be modified by collective effects. We use a 
hard thermal loop effective gluon propagator. The gluon propagator is diagonal in color, and color indices 

will be omitted. We work in the strict Coulomb gauge and decompose the gluon propagator into transverse 
and longitudinal pieces: 

-D% /A (L) = Pf(L) A T R/A (L) + P£»(L) A L R/A (L) (2.3) 
The projectors are defined as [231 : 

w)=(S V+w) (2 - 4a) 

L 2 f I 



TW = T o o ) ( 2 - 4b ) 
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and coefficients are: 



A T / L (L) = 



L2-n T/L (L) 



(2.5) 



The quantities ILjy^L) are defined by choosing the decomposition of the gluon polarization tensor that 
satisfies the usual Schwinger-Dyson equation: 



n""(L) - Pf(L)Il T (L) + Pf(L)IL T (L) = ^((^(L))- 1 - (D^(L))~ l ) 



(2.6) 



where Dq V (L) is obtained from H2.3J1 by setting the polarization tensor to zero. We note that l|2. 6J1 holds 
separately for retarded and advanced components. The subscripts R/A have been omitted to simplify the 
notation. We use the 1-loop HTL results 0] HEj : 



Ht(Ix, l) = 3 m 



x 

T 



x(l — x 2 ) fx + 1 
In 1 



U L {lx,l) = 3m 2 (l-x 2 ) 



x , 

1 In 

2 



x-1 
x + 1 
x- 1 



(2.7a) 
(2.7b) 



flavors. The spectral function is defined as 



where x = Iq/1 and = g 2 T 2 [N + Nf /2]/9 is the soft gluon thermal mass in a SU(N) gauge theory with Nf 



PT/L (loJ) = A T B /L (L)-A\^(L). 



T/L, 



(2.8) 



3 Calculation 

We want to calculate the trace of the imaginary part of the retarded photon polarization tensor which is 
given by 

Im(IT^ = | ((Hfl)^ - (ITOg 

The diagrams to be evaluated are shown in Fig. [1]. Since we consider the production of a static external 
photon (Q M = (^0)0)), with qq > 0, all the azimuthal angular integrations can be carried out trivially and 
the integration measure is rewritten as, 

d 4 P [ d 4 L = 8tt 2 / dl [ dll 2 [ dy [ d Po [ dp p 2 (3.1) 



where y = I ■ p is related to the angle between the gluon and quark momenta. We simplify the notation by 
defining the momentum variable R = P + Q. 

3.1 Vertex diagram 

The vertex contribution (Fig. [la]) to the imaginary part of the photon polarization tensor can be written 

as (2E1EZ1EH1: 



2 » 2 r d i P r d ± L 



X 



X 



D R (P)D R (P + L)- Da{P)D a {P + L) 

D R (R)D R (R + L) - D A {R)D A (R + L)] (3.2) 



where we have kept only terms proportional to the gluon thermal factor ub(Io) since we are interested in a 
possible enhancement from the soft part of the gluon phase space. Rewriting this expression using H2.2J) we 
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obtain, 



imn 



(vertex) " 
R u 



N c C F e 2 g 2 



d*_P_ 

(2^ 



[n F (r ) - n F (po)] p(P) 



d 4 L 
(2tt)4 



Tl (verte X ) nB y o)nP o 



x V(L + P) p(R + L)V(R)+p(R)V(R + L) 

where we have defined 

■R^= PT (L)Pf(L)+p L (L)P^(L), 
and the Dirac trace is calculated to be: 



(3.3) 



Tr 



(vertex) 

pa — 



Tr [f> r Hl P (H + 1) lp (f + t) la] (3.4) 
-8(2£ -PL-R + L- RP 2 + L ■ PR 2 + P 2 R 2 - L 2 R ■ P)g pa 
+16(R 2 (L P + P p )P a + P 2 {L P + R p )R a + L ■ RP p {L a + P a - R a ) 
+L ■ PR p {L a +R a - P a ) - L 2 P p R a - R ■ P{P p R a + (L p + P p ){L a + R a )) . 

By analogy with the zero temperature calculation, the first term in eq. (|3.3(l is called the real cut (the fermion 
propagators P and R + L are cut) while the second term is called the virtual cut (P and R are cut). There 
are two additional cuts: one where the fermion propagators carrying momenta P + L are R are cut, and 
one where the propagators carrying momenta P + L and R + L are cut. Using symmetry arguments one 
can show that these two terms contribute an overall factor of two which has been included in (|3.3|) . The 
contraction P^ L r Ti p v ^ rtex ^ is straightforward to perform. We look at the longitudinal and transverse parts 
separately. We obtain, 

p pa Tl ,(verte x) = g(1 _ ^ [(fo2(2p2 + l{pQX + py) _ q Q l(-2p 2 X + p Q (l(l - x 2 ) + 2py)) - 2l 2 p 2 (l - y 2 )] 

pP^verte.) = ^ { q 2 {p 2 {1 _ y2) + ^ _ py)) + + ^) _ 2 pxy) ~ 2lp(p y - px)) 

- l 2 p{4 Po xy - p(l + x 2 ){\ + y 2 )) } (3.5) 

where we have used x :— and set P 2 = in anticipation of the fact that the factor p(P) in (|3.3|l will 
give P 2 = when 1|2.1|) is used. 



3.2 Self-energy graph 

The expression for the self-energy diagrams is similar. We obtain the following expression for the diagram 
of Fig. [lb]: 



im n 



self(b)P 
R u. 



~N C C, 



e 2 g 2 



d A P 



[n F (r ) - n F (p )} 



d*L 



Tr p s f nb)) n B (l )KP° 



x D R (P) - D A (P)\ \{D R {R)) 2 D R {R + L)-{D A {R)) 2 D A {R + L)\ (3.6) 
where, as before, we have kept only the terms proportional to hb((o). Rewriting using (|2.2|l we obtain, 



Im n 



R u 



N c C F e 2 g 2 



d^P 
(2^ 

x 



[n F (r ) - n F (p )] p{P) 



d i L 



(2tt) 4 pa 



V(R + L) V(R) p(R) + (P(R)) 2 -p(R + L) 



(3.7) 



As in the case of the vertex diagram, the first term in the square bracket corresponds to the real cut, and 
the second term corresponds to the virtual one. The Dirac trace is: 



rp self(b) 
- L1 pa 



= Tr[f ^ Hip (fi + £)iafilp] 

= 8(g pa (R 2 {R ■ P — L ■ P) + 2L ■ RR ■ P) + 2{L p + R p ){R 2 P a - 2R ■ PR a )) 



(3.8) 
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The expression for the second self-energy diagram (Fig. \3p]) can be obtained in the same way and it is taken 
into account by multiplying eq. 1)3. 7|) by a factor two. Contracting P^ L Tr s p ^^ we obtain, 

pP^sejm = 8(1 _ x 2 )q 2 {qQpo + l{pQX _ py)) 

p^T^m = -Wql(p oqo + (l(p x+py)+p 2 (l + y 2 )) (3.9) 

where we have set P 2 = as before. 

The first step in the calculation of (|3 . 3|) and i|3.7[l is to perform the po integration using (|2.1|) . The 
choice qo > means that only n = 1 is kinematically allowed. The second step is to use the delta functions 
contained in the spectral functions p(R + L) and p(R) to perform the p integrations in the real and virtual 
terms respectively. The structure of the real cut is the same for both vertex and self energy diagrams. 
Making use of (|2.1|l we have, 

in , t\ n xi \ ik + qo) 2 - 1 2 , Q 1fV . 

p(R + L) = — — -6(p - p rea i) ; p roa i = — — — — (3.10) 

lo + qo + ly 2(l + qo + ly) 

Similarly, we attack the virtual cut by using l|2.1|l to obtain (remember that r = q a + Po and r = p in static 
photon kinematics), 

P(R) = —5{p - Pvirtual) ; ^virtual = 7T 

qo 2 

This expression can be used directly to do the p integral in the virtual part of the vertex diagram. The 
structure of the virtual part of the self energy diagram is more complicated: there is a singularity because of 
the presence of the factor 7 , (i?) p{R). This singularity can be regulated in the standard way We write 



I \ _ d ( 1 

dx 



(3.12) 



Taking the imaginary part of this equation gives, 



Using this expression we obtain, 

d_ 

H dp 



V{R) p(R) = — —6(p - ^virtual) (3.14) 



and perform an integration by parts to remove the derivative from the delta function. 

Once the variables po and p have been integrated over, we can simplify the result by performing an 
expansion in Zo/<7o an d l/qo- However, as shown below, one finds that the leading order term cancels between 
the real and virtual contributions to each diagram, and thus we must expand to next-to-leading order. The 
fact that p r eai is a function of Iq, I and qo means that the contributions from the real cuts will involve an 
expansion of the (npiro) — np(po)) thermal factors. We use the notation: 

rip = e p qo j 2 T i > n F = Q; n' F = n F (l - n F )(l - 2n F ) (3.15) 

After the expansion has been performed, the resulting expression can be integrated over y by hand. We look 
at longitudinal and transverse pieces separately. 

4 Preliminary Results 

We list below the results for the trace of the imaginary part of the retarded polarization tensor. To simplify 
the notation, a common factor has been extracted from each term: 

¥kC = -^£ l N c C F J dll J dxn B (xl)p L/T (lx,l) (4.1) 
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In addition, we have used the notation: 

x + 1 



LOG = In 



x-l 



X = — - 



We obtain: 



VertexX irtual - (1 - 2n F ) — — (g 2 LOG - 2l 2 {2x + (1 - x 2 )LOG) 
Vertex^' = — (1 — 2n F ) ~~~~~~ [?oLOG — 2l 2 (x + (1 — x 2 )LOG] 

Vertex L = -l 2 (l - x 2 ) ^2(1 - 2n F ) + - jcLOGJn'^ 



Self - EnergyX irtual = 2 ( 1 - 2n F ) q (q + 2lx) 



Self - Energy r L cal = -2(1 - 2n F ) [q 2 + 2lq x - l 2 (l - x 2 )] - -Z 2 ||(l - x 2 )(l + 
Sclf-Encrgy L = / 2 (1 - x 2 ) ( 2 (1 - 2n F ) \ ||(1 + 



LONG = -Ip||(l-x 2 )(2- a ;LOG + — X 
2 -t qo 



Vcrtex T rtual = (1 - 2n F ) - (q 2 - 2l 2 )(2x + (1 - a; 2 )LOG) 

Vertex™" 1 = -(1 - 2n F ) - [q 2 {2x + (1 - a; 2 )LOG) - 2Z 2 (1 - a; 2 )LOG] 



-^l 2 Q(1- 32; 2 )-^(1-. 2 )LOg) n' F 
VertexT = -I 2 (a(1 - 2n F ) + || Q(l - 3x 2 ) - ^(1 - x 2 )LOG^ n'^j 



Self - Energy^* 1 " 2 ' = -2(1 - 2n F )[ g 2 (2 - xLOG) + l 2 x(2x - (1 + x 2 )LOG] 
Self - Energy^ = 2(1 - 2n F ) [(2 - xLOG)(g 2 + l 2 {\ + x 2 ))} + ^||; 2 n F 



Soir-Enor-v, = r [A{l-2n F )+ l -^n' F 



TRANS = - ,2 ||x(2x + (l-x 2 )LOG)n F 
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Before proceeding, we make the following observations about these results: 

[1]: In each of the four cases: (vertex/self-energy) x (longitudinal/transverse), the leading order term in the 
l/qo expansion cancels between real and virtual contributions and the final result is of order I 2 . 
[2]: For both the longitudinal and transverse pieces the term proportional to (1 — 2uf) cancels between the 
self energy and vertex contributions. 

[3]: The longitudinal propagator has the form Al(L) ~ 1/(1 — a; 2 ), as can be seen from (|2.5[) and 1)2.70 . 
Consequently, the terms (|4.5|l and 14.611 exhibit co-linear singularities in the form of a factor 1/(1 — x 2 ) which 
diverges when x — > ±1. When the two expressions are combined however, the sum i|4.7|) is well behaved. 

As a check on these results, we sketch in the appendix an alternative derivation of the transverse gluon 
contribution. 



5 Iq and I integrals 

In this section we perform the two remaining integrals over Iq and I. From 1)4. 1|) . (|4.8)l and (|4.15)l and using 
the approximation hb{Ix) ~ T/(lx) we have, 



LONG = 



(27T)' 



-N c C F n F I L 



T 



TRANS = ~ 77 ^rjN c C F % n' F I T 



(27T) 



where, 



fOO J poo 

Il= — dll 2 p L (lx,l) 

J-oo x JO 
r°° rj T f°° 

i T = ££ / dii 2 PT (ix,i) 

J-oo x Jo 



T 



(l-* 2 )> 



x + 1 



a: - 1 



x 2 + (l-x 2 )-ln 



- 1 
x + 1 



a - 1 



(5.1) 
(5.2) 



It is interesting to remark that the terms in square brackets in the integrands are proportional to Relli and 
RellT respectively. 



5.1 Longitudinal Integral 

We consider first the longitudinal integral (|5.1|) . We divide the integral into two pieces, one with and one 
without the logarithm: 



n 



dx 



dll 2 p L {lx,l){\- x 2 ) 



-co * JO 

CO 



dx 



dll 1 p L (lx,l) (1 - x 2 ) - In 
x ,/ 2 



x + 1 



x-1 



The calculation of the first piece is standard |2H]. One obtains 

If = %/37r 2 m 



(5.3) 
(5.4) 

(5.5) 



To calculate l\ we use the method proposed in [23 • We use an integral representation for the logarithm and 
write, 



In 



x + 1 



x- 1 



dx'V 



1 



(5.6) 



where the principal part is defined as 



1 f 1 



2 V x + ie x — it 
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The first step is to perform the integral 

dx±(l~x 2 ) p L (l x J) V (—L^j (5.7) 

via the standard dispersion relation (see 29 ) keeping track of the contribution of the large circle at infinity. 
We obtain 



l\ = -TrRe / dy j d'l [ I 



-1 Jo 

and performing the integral over / we have, 



R r 1 

I b = ^-7rm s y da;'ReVQi(a;') (ln(VQi(x')) - ln(-VQi(^)) 



(5.9) 



Cutting the complex plane along the negative real axis for both the square root function and the logarithm 
we find, 

(ln(VOiF)) - H-^Qi(x'))) = -iTtSgn(x') 

which gives, 

/& = ^-vr 2 m 3 Im y dx'sgn(a; , )v / Qi(2 ; ') (5-10) 

We do the x' integral by writing f(z) = \J Qi{z) and choosing the contour shown in Fig. (J2J). We obtain, 

h = -Tr 3 m g /2 (5.11) 




Figure 2: Contour in the complex x' plane 
Combining (|5.5|) and (|5.11|) we have, 

/ L =7r 2 m 3 (V3-|) (5.12) 

5.2 Transverse Integral 

In this section we consider the transverse integral l|5.2() . The £- integral can be performed in a way similar to 
the one above, but the result exhibits an ultraviolet divergence in the Z-integral. To understand this result, 
we separate the gluon spectral function into the pole piece and the cut piece (see [291). 

The pole piece carries the ultraviolet divergence. This divergence is not physical: it is a consequence 
of performing an infrared expansion on an integral that is dominated by the ultraviolet part of the phase 
space (the expansion in l/qo used above is not valid in this case). But the pole part is just the contribution 
calculated in 9 (up to higher order terms since the transverse gluon has a small thermal mass), the result 
of which is given in eqs. I|1.3|) . i|1.4fl . 
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The cut part of the spectral function does not produce a divergence. We have, 



r 1 At r°° 
t= / ax / apfo, * 

J-i x Jo 



x 2 + (l-x 2 )^-ln 



1 



(5.13) 



where 



0r{x,t) 



-2 Im n T 



(l 2 {x 2 - 1) - Re n T ) 2 + (Im n T ) 2 



and ILt is given in l|2.7a|l . Since the integration range does not extend from positive to negative infinity we 
cannot use a sum-rule-like trick to perform the integral. We proceed as follows. The J dl integral can be 
done by hand. The only mass scale in the integrand is m g and a simple dimensional analysis shows that the 
J dx must be of order m g ~ gT. The remaining integral is done numerically to obtain, 



rCUt 



-6.67 m 



(5.14) 



which is of the same order as the longitudinal contribution. 



6 Conclusions 

In this note we have computed the lepton pair production rate by a static heavy photon at two-loop or- 
der (Fig. PP) and focused on the infrared region for the exchanged gluon. The result is free of infrared 
and collinear divergences, in agreement with a previous study in the case where the external particle was a 
scalar As in 1321. we find the contribution of the longitudinal gluon to be of order e 2 g 3 . However, our 
results differ for the transverse sector. The authors of |22] find that the contribution of the transverse gluon 
is of order e 2 g A . Our results are as follows: 

[1] The infrared approximation used in this computation is valid for the "cut" part {x 2 < 1) of the gluon 
spectral function. The contribution for "cut" part is of order e 2 g 3 . 

[2] The infrared region of the gluon phase space does not give the dominant contribution to the pole piece of 
the gluon spectral function. The pole piece is dominated by large momenta. This means that thermal effects 
can be neglected and we can use the result obtained with a bare gluon propagator which was calculated 
previously UJ. 

The physical picture of the calculated contributions is the following. The one-loop contribution is the 
annihilation of a quark-antiquark pair, each parton carrying half of the energy of the static photon, namely 
qo/2. The dominant two-loop contribution describes the scattering of one of the annihilating quark in the 
medium with an exchange of small momentum of order I ~ m g ~ gT. Interestingly, these processes are 
similar to those which dominate the production of the small mass, high momentum (of 0(T) or higher) 
lepton pairs, i.e. the other extreme kinematical domain |3()| . In the static case, higher order loop diagrams 
may modify the result, as was found to be the case for energetic photons with a small or vanishing invariant 
mass |3"Tll3"2] . 

We summarize below the result for Imll/j^ up to order g 3 in perturbation theory : 



ten R » = aN c <fi[l-2n F ] 



1.62 UF (l-n F )^f 



(6.1) 



The complete result <|6.1[1 shows that the processes considered in this paper can be relatively large when 
qo/T is not too large. This suggests that these processes may be relevant for lepton pair production from 
static heavy photons at temperatures where we expect a quark-gluon plasma to be formed |33l 1341 135| . 



7 Appendix 

As a check of the above results we have derived the transverse gluon contribution using an alternative method. 
The J dpo integration is done using Cauchy's theorem in the complex po plane. The J dy integration is 
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performed second, and then the imaginary part is taken. The J dp integration is done last. In addition, 
vertex and self-energy matrix elements are combined before carrying out the integral so that a partial 
cancellation of terms occurs before integration. We start from 1)3. 2|) and (|3 .(j|> . In the po complex plane, the 
integrand has poles above and below the real axis from the propagators, as well as poles on the imaginary 
axis from the statistical factors n F (ro) and n F (pa). First we note that when we close the contour in the 
upper-half or the lower-half po plane, the contribution from the poles on the imaginary axis drop out because 
the difference between the retarded and the advanced propagators vanishes. Second, we use the fact that 
terms containing only retarded or advanced propagators do not contribute. The product of propagators can 
thus be simplified so that we obtain: 



Im n 



self(b)» 
R u 



N c C F e 2 g 2 Re 



d 4 P 



[n F (r ) - n F {p )} 



d 4 L 



x D A (P)(D R (R)) 2 D R {R + L) 



Tr(f^ b »n B (Zo)^ po 



(7.1) 



Imll 



(vertex)" 

R u ~ 



N c C F e 2 g 2 Re 



d 4 P 



[n F (r ) - n F (p )} 



d^L 



(2tt) 4 l v u/ -yrvnj 
D R (P)D R (P + L)D A {R)D A (R + L 



Tr<£ ertex) n B (l )K p ° 



(7.2) 



In the expressions above it is understood that the poles in the statistical weights should be ignored when 
doing the po integration in the complex plane. Only the poles from the propagators contribute. 

We can now add vertex and self-energy contributions so that the full expresion of the trace on the quark 
loop multiplied by the appropriate propagators is: 

Tr pcr ■ propagators = 



2L 2 (RPR a + PPP a ) - AQ 2 RPP a + gP a L 2 (Q 2 + L 2 - P 2 - R 2 - (P + L) 2 - (R + Lf) 



P 2 R 2 A (P + L) 2 R (R + L) 2 



Q 2 



-4PPP a + g pr7 {(P + L) 2 - L 2 

{pWrW+Wr 

1 1 \ ( 1 



-4RPR a + g pa {(R + L) 2 - L 2 



P 2 R (R 2 A ) 2 (R 



p2 
1 R 



{P + L)] 



R 2 A 



(R + L)\ 



fJ 



per 



L? A 
2QL 



(7.3) 



2QL 



P 2 R 2 A (P + L) 2 R P 2 R A (R + L) 2 A 



where we have used the notation P R 



(po+ief—p 2 , etc. As discussed in (see Eq.(3.12) of this reference), 
this expression already exhibits some cancellations due to gauge invariance. Further simplifications are 
possible in our case since we are working in the limit of large q 2 = Q 2 and small {Iq, I 2 , L 2 }. At this point 
we consider only transverse gluons (longitudinal gluons can be treated using the same method) and identify 
the terms that contribute to (|7.3(l to leading order. It has been checked explicitly that all other terms are 
subdominant. We need to calculate: 



Ti = e 2 g 2 N c C F Re 
T 2 = e 2 g 2 N c C F Re 



d 4 P 



(2^)4 
d 4 P 

W) 



[n F (r ) - n F (p )} 
4 [n F (r ) - n F (p )} 



d*L 



(2tt)4 
d 4 L 



n B (lo)pT{L) 



(27T) 



4 n B (lo)pT(L) p URl)2{R 



16 Q 2 RPP a Pf(L) 
P 2 R R A (P + L) 2 R (R + L) A 
16 Q 2 RPR a P pa {L) 



L? A 



(7.4) 



and an equation similar to T 2 with P and R interchanged. We discuss below the strategy for calculating the 
necessary integrals. Consider, for example, a term of the form: 



Int = Re / p dp dy dpo 



P 2 R (P + L) 2 R R 2 A (R + L)< 



(7.5) 



where the integrals over l and I have been supressed (see H3.ll) ). These integrals must be calculated following 
the method of Sec. and are not discussed in this appendix. We close the p contour in the lower half-plane. 
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There are contributions from the poles of the factors Pj| and (P + L) 2 R . In the small {Iq, 1} limit the 
contributions from the poles of the factor (P + L) 2 R will have the same form as the contribution from the 
poles of Pft. Thus we will consider only the poles po — ±p — ie from the factor Pjj, and multiply the result 
by a factor of two. This procedure gives rise to an integral of the form, 

Re [ dpo E±Pf> = VRef-P(np-z £ ) 



P R t± V P 

= --TmF(-p-ie), (7.6) 
P 

where we have neglected, in the last line, the pole which will not lead to an imaginary part in the function 
F. Using (|7.6[1 we obtain, 

* I 1 1 (77) 

R\ q (qo-2p-iey (P + L)\ c-by' (R + L) 2 A a - by 

where 

a - b = (q + l Q + I) (q - 2p + l - I) - ie, a + b= (q + l - I) (q - 2p + l + I) - ie 

c-b= (l + l) (-2p + l -l), c + b= (l -l) (-2p + l + l). (7.8) 

We note that when the retarded propagator 1/(P + L) 2 R is evaluated at the pole of the retarded function 
1/Pr, the ie piece cancels exactly and thus the variables b — c and b + c in (|7.8I) do not contain an ie 
prescription. Carrying out the angular integration on y, one obtains: 

x dV (a- by) (c - by) ~ 2q a pl (q - 2p + 2l - ie) (o + ft) ~ ^ (j+l) ) ' (7 ' 9) 
Combining lfT5|l . l(7~B|) . ifTTjl and O we have, 

Int (7.10) 
= — 2tt / dpplm 



where 



Zqopl Qa(qo -2p- ie) (q - 2p + 2l - ie) 



1 1 / f 1 \ ( 1 



(qo-2p-ie)(q Q -2p+2l -ie) 2l \ \q -2pj \q -2p + 2l a 



ITT 



+ — (S(qo-2p)-S(q a -2p + 2l )). (7.11) 

There is also an imaginary piece to the first logarithm in 17.100 . We use the relation ln(x — ie) — hi(|ir|) — 
iirQ(—x) to separate real and imaginary pieces and obtain, 

\a + bj J \ \q + l -lj \q a -2p + l Q + lJ J J( qo+ i -i/2 

To complete the evaluation of l|7.10|l we take the product of H7.11|l and (|7.12|) and extract the imaginary part. 
There are two contributions: one contains the imaginary part of (|7.11|) which contains a delta function that 
allows us to do the J p integral, and the second contains the imaginary part of l|7.12|l . This factor contains 
an integral in p phase-space that is restricted to a small interval around qo/2 which can be calculated by 
expanding around the value p = qo/2. This procedure leads to the appearance of a factor n" F as defined in 
(|3.15|) . The J dp integration can then easily be done. The calculation of terms like T 2 is carried out in the 
same way, the double pole being handled via an integration by parts on the momentum variable after taking 
the residue of the single pole via the contour integration. Finally we obtain: 



Ti = J^y[ N c C Pf n F J — I dl I z Pt(x,1) 



x 2 -l + (l-x 2 )^-\n 



(7.13) 
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T 2 = ±^ C,!„£ /f |<fl lVr(*,0j- (7-14) 
Adding l|7.13[) and twice l|7.14|) gives the full transverse contribution as calculated in l|5.2l) . 
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